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Abstract 

It is stressed that a measurement of the electric dipole amplitude for direct photon emis- 
sion m — > 7r ± 7r 7 decays through its interference with inner bremsstrahlung is important 
for differentiating among various models. Effects of amplitude CP violation in the radiative 
decays of the charged kaon are analyzed in the Standard Model in conjunction with the 
large N c approach. We point out that gluon and electromagnetic penguin contributions to 
the CP- violating asymmetry between the Dalitz plots of K ± — > 7r ± 7r°7 are of equal weight. 
The magnitude of CP asymmetry ranges from 2 x 1CU 6 to 1 x 10~ 5 when the photon energy 
in the kaon rest frame varies from 50 MeV to 170 MeV. 
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In a recent paper [1] we have studied CP violation in the radiative kaon decay Kl — > 
7T + 7r~7. We conclude that the direct CP-violating effect originating from the electromagnetic 
penguin diagram is only of order (1CT 3 — lCT 4 )e, depending on the region of the Dalitz 
plot under consideration. On the contrary, it has been advocated that direct CP-violating 
asymmetry in the radiative decays K ± — > 7r ± 7r°7 defined by 

T(K+ -> 7r+7r°7) - T(K~ -> 7r-7r 7) 

r ~ T(K+ -> 7T + 7T 7) + r(^- -> 7T-7T 7) ' ^ ^ 

arising from the same electromagnetic penguin mechanism can be large enough for experi- 
mental interest; explicitly, A r < 9 x 10~ 4 is obtained in Ref.[2]. If this estimate is correct, 
it will be on the verge of the capability of the factory DA$NE [3] , and could be detected 
at future high-statistics facilities. The purpose of this short Letter is to re-examine this 
CP- violating effect in the Standard Model in conjunction with the l/N c approach. 

The general amplitude of the decay K + (k) — > -K + (p + )-K°(p)'y(q, e) is of the form 

-> tt+A) = -eA(K + -> tt+tt ) (?±-^ - e< (2) 

\p+-q p-q J 

+M[iee llvpa p> x + p»q»e°\e^ + Ee[(p + ■ e)(p ■ q) - {p ■ e)(p+ ■ q)]e^, 

where we have included the isospin phase shifts Sj with / being the isospin of the two pions 
and J the total angular momentum, which are necessary for generating the decay-rate CP 
asymmetry. To the leading multipole expansion, the direct emission (DE) amplitude M 
corresponds to a magnetic dipole (Ml) transition, while E an electric dipole (El) transition. 
From time to time this decay mode has received a constant attention both theoretically [2-14] 
and experimentally [15-17]. Outside of the framework of chiral perturbation theory (ChPT), 
various techniques such as the short-distance effective weak Hamiltonian, the vector-meson- 
dominance model and current algebra have been employed to study the DE of K + — > 7r + 7r°7 
decay. These methods are plagued with some fundamental problems. For example, the 
short- distance effective Hamiltonian utilized in Refs.[5,7] does not explicitly couple to the 
external photon field. Hence, after the usage of the factorization approximation, one has to 
appeal to the soft-pion theorem to evaluate the matrix element (7r°7|s7^(l — r y^)u\K + ) ) for 
instance. However, it has been shown [8] that the soft-pion technique is no longer applicable 
to the magnetic transition amplitude as in the case of n° — > 77. Also, it is known that the 
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short- distance effective weak Hamiltonian is far from being adequate to describe the kaon 
AI = \ rule. 

In ChPT, the most general p CP-invariant AS* = 1 non-anomalous electroweak chiral 
Lagrangian with one external photon field which satisfies the constraints of chiral and CPS 
symmetry has the expression [18] 
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Cf-Lom = i^g^[u{l\{X^L v Q)-ru 2 Tr{X,L v QL^ (3) 
+ ousTriXeUR^RvQtf) + u 4 Ti {X^UQR^RM^) + uj^t(X 6 UR v QR^) 



for normal intrinsic parity transitions, while anomalous Lagrangian terms for the odd intrin- 
sic parity sector are [10,12] 

' 2 ' 
2 ' 

T 2 

, J IT 



^ = ia[^)g 8 eF^Tr(QL,)Tr(X 6 L v ) 



+ ib{-\ g 8 eF^Tr(QRjTr(X 6 L v ) 



+ ic(j^J g 8 eF^Tr(x 6 [UQU\ L„L„]) , (4) 

where F^ = e^F^, Q = diag(|, -§, -§), L M = (D^U)U^ with D^U = d^U-ieA^Q, U], 
Rf, = U^D^U), and U = exp (2i£) with f w = 132 MeV, = ^ a X a . In Eqs.(3) and (4), 
gs is the octet weak coupling constant appearing in the lowest order CP-invariant AS = 1 
weak chiral Lagrangian CyJ = —g 8 Ti{X & L^L^) and is fixed to be [19] 

# 8 = -0.26 x 10~ 8 m^ (5) 

from the experimental measurement of K° — > nn decay rates. The coupling constants 
Ui, a, b, c depend on the choice of the renormalization scale ji as divergences of chiral 
loops are absorbed by the counterterms which have the same structure as £^n-anom an d 
^anom 1 - Therefore, those coupling constants in principle can be determined only empirically 
from various low-energy hadronic processes. However, in the limit of large iV c (N c being 
the number of quark color degrees of freedom), these couplings become \i independent and 
are theoretically manageable at least to the zeroth order of a s [10]. It is found that in the 
large- N c approach [10] 

N c N c 
UJ i = UJ i = 7^> uj 3 = u 4 = uj 5 = 0, a = 2b = 4c= (6) 



Note that the couplings a, b and c are determined by chiral anomalies and hence are free of 
gluonic corrections in the large N c limit. 

There are two different contributions to the direct emission of K ± — > it^-n ^: contact- 
term contributions induced by £^n-anom and C^^, and three long-distance pole diagrams 
with the 7T7T7r7 vertex governed by the anomalous Wess-Zumino-Witten term. As shown in 
Ref.[10], the results are (Since we are working in the leading order in 1/N C expansion, chiral 
loops can be neglected.) 

EiK* -> 7T ± A) = Contact = 4f (2), 



M (K ± - n ± 7r° 1 ) = M contact + M pole = T^§(2 + 3). (7) 



The constructive interference between pole and direct-transition Ml amplitudes for K ± — > 
7r ± 7r°7 decays is a prominent feature different from the decay Kl — > 7r + 7r~7 where a large 
and destructive interference for Ml transitions is required to explain the data [1,10]. Ex- 
perimentally, the DE rates are extracted in the charged-pion kinetic energy range of 55 to 
90 MeV [15-17]. With this experimental condition, the branching ratio of direct emission is 
given by [10] 

BriK* -> 7r ± 7r°7) D E = 1.32 x 10 5 (|£| 2 + |M| 2 ) GeV 6 

= 2.02 x 10" 5 , (8) 

which is in agreement with the experimental values 

' (1.56 ± 0.35 ± 0.5) x 10- 5 [1972] (Ref.[15]), 
BriK* -> 7T ± 7r°7) DE = | (2.3 ± 3.2) x 10" 5 [1976] (Ref.[16]), (9) 

. (2.05 ± 0.46±g;ii) x 10" 5 [1987] (Ref.[17]). 

Previous calculations [5,7,8] based on the short- distance effective weak Hamiltonian predict 
a smaller branching ratio. This is attributed to the fact that, as noted in passing, only short- 
distance corrections to the Wilson coefficient functions are taken into account in the approach 
of the effective weak Hamiltonian, which are not sufficient to explain the AI = | rule in 
kaon decays (see Ref.[19] for a review). We note that although the DE rate is dominated by 
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magnetic transitions [] due to additional constructive contributions from the pole diagrams, 
the El contribution is nevertheless not negligible. Experimentally, the DE electric dipole 
amplitude can be measured from the interference of inner bremsstrahlung (IB) with El 
transitions. Thus far, there is only one experiment (done 2 decades ago) measuring this 
interference and only a limit is obtained [15]. Evidently, a measurement of the El amplitude 
is important for understanding the underlying mechanism of K — > 7r ± 7T°7 decays and their 
direct CP violation. 

We next turn to examine CP-violating effects in the decays if — > 7r ± 7r°7. The simplest 
way of observing CP nonconservation is through the measurement of the decay-rate CP 
asymmetry parameter Ar as defined in Eq.(l). If photon polarization is not measured, there 
is no interference between El and Ml amplitudes. Consequently, when photon polarizations 
are summed over, a nonvanishing Ar must arise from the interference of IB with the El 
amplitude of DE. Although CP-violating asymmetry has been studied intensively in late 
1960's [20], a modern analysis in the framework of the Standard Model was carried out only 
recently in Refs.[2] and [7]. However, the result is somewhat controversial: While Ar arising 
from the gluon penguin diagram is estimated to be of order 10~ 6 by McGuigan and Sanda 
[7], it is claimed by Dib and Peccei [2] that a large CP asymmetry can be induced from 
the electromagnetic penguin diagram, namely Ar < 9 x 1CT 4 . Naively, it is expected that 
amplitude CP violation coming from the electromagnetic penguin diagram with a photon 
radiated from the loop quark or from the W boson is of equal weight as that from the QCD 
penguin diagram with a photon emitted from the external quark lines or from the W boson. 
Therefore, a resolution of this discrepancy is called for. 

We will follow Ref. [7] to consider CP asymmetry between the Dalitz plots of K + — > 7r + 7r°7 
1 Precisely, \E/M\ 2 = 0.16 is predicted in the large-iV c approach. Most earlier calculations yield even 
smaller ratio for \E/M\ 2 . In Ref. [14] this ratio is calculated to be 5.1 x 10 -3 . In ChPT, the coupling constants 
of Anm-anoni an d ^tmym are expected to be of the same order of magnitude [see Eq.(6)]. Therefore, unlike 
the decay Ki — > 7r + 7r~7, it seems to us that there is no reason to have a severe suppression on the electric 
dipole amplitude of DE in K — ► 7r ± 7r°7 decays. 
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and K — > 7r 7r°7 rather than in the total decay rates 

. \A(K+ ^7T+7r 7 )| 2 - \A(K~ ^7T-7r 7 )| 2 



(10) 



-> 7r+7r°7)| 2 + -> 7r-7r 7)| 2 ' 

so that a larger CP asymmetry can be obtained in certain particular regions of the Dalitz 
plot. Since under CPT invariance 

A(K~ -> tt-ttV) = -eA(X- -> tTtt ) - ^] e< (11) 

\p--q p-qj 

-M*[iee^ P tp v q p e°\e i5 i + £*e[(P- • e)(p • q) - (p • e)(p_ • g)]e^, 

it follows that 

a = — 2|^7|^sin0 jE ?siii(^ - jg) , . 

A(K+^7r+7r o )/m^ + 2|E|zcos0 E cos(5 1 1 -5 2 ) + (|E| 2 + |M| 2 )2 2 ' 1 J 

where = — \E\e l ^ E [recall that our E is negative; see Eqs.(5) and (7)] and z = (q • p±)(q ■ 

p)/m 4 K . 

The main task is to estimate the CP-odd phase 4>e of the El amplitude. There are two 
different contributions to the imaginary part of E: one from the gluon penguin diagram, 
and the other from the electromagnetic penguin diagram. As to the former, following the 
prescription presented in Ref. [1], we find in the l/N c approach that 

where g' 8 is the CP-violating coupling constant appearing in the lowest order CP-odd AS* = 1 
weak chiral Lagrangian £^ = —ig' 8 Tr(\ 7 L IJj L fJ ') and is dominated by the short-distance QCD 
penguin diagram. Hence, 

(sm0£) gluon = — = ^-r-, (14) 

where A = A(K° — > itit(I = 0)), and use of 

A(K 2 -» = 0)) 

A(K, - n7r(I = 0)) <? 8 1 J 

has been made. The calculation of Iim4 in the Standard Model is standard and is given by 
[21] 

lmA = -^(ImA t )y 6 ((7r7r) /=0 |g 6 |ir°), (16) 
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E = zhiHs + tfs), (13) 



where y% = Imce/Imr, r = —\ t /\ u , A« = V* s Vid, Qe is a penguin operator and cq is the 
corresponding Wilson coefficient. The K — tttt matrix element of Qq evaluated in the large 
N c approach is known to be [19] 

((7T7r) /= o|g 6 |ir°) = -i±V3Uv 2 m K-< ^ (17) 

X 

where 

m l _ m K+ _ m \« , 18) 



m u + m d m u + m s m d + m s 

characterizes the quark order parameter (qq) , and A x ~ 1 GeV is a chiral symmetry breaking 
scale. Since ReA = i4.69 x 10~ 7 GeV [19] and 

lm(V t * s V td ) ^ si 3 s 23 sin 5 13 (19) 

in the Chau-Keung parametrization of the quark mixing mixing matrix [22], we find numer- 
ically 

(sin B ) g i uon = -8.1 x Kn 5 sin<5i3, (20) 

where uses have been made of m s = 175 MeV, s 23 = 0.044, S13/S23 = 0.1, and y 6 = —0.057 
for m t = 150 GeV [21]. 

Following Ref.[l], the DE amplitude of K + — > 7r + 7r°7 induced by the electromagnetic 
penguin diagram is given by 

A(K L -> Tr+Tr-7)^ = iG F T ^Im(V; s V td )F(x t )(7t + 7r° 1 \Q T \K + } } (21) 



with 



Q T = i[m s sa f j lU (l-"f5)d + m d sa tJ , u (l + j5)d]F^, 
(8x 2 + 5x - 7)x (3s - 2)x 2 

r (X) = ; r- ; ill X , 22 

v ; 12 (x-1 3 2 (x-1 4 V ; 



2 The experimental values of K — > tttt amplitudes are usually expressed in terms of real numbers. However, 
model calculations show that the amplitude of K — ► 2ir contains a factor of i. This means that Re^lo given 
in Ref.[19] should be multiplied by a factor of i when compared with (tttt\Qq\K q ). 
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and Xt = m\/Myy- By working out the chiral realization of the tensor operator Qt as in 
Ref. [2] (see also Ref.[l]), we obtain 

Eem = i o ~S~T7? F ^ ( g i3^23 sin S 13 ). (23) 
It follows from Eqs.(23) and (7) that 

(sm(f) E ) em = ^f em = ^ GFms ^ F(x t )(s 13 s 23 sin 5 13 ). (24) 
E 8 gs 

Numerically, 

(sin <j) E ) em = -6.0 x KT 5 sin 5 13 , (25) 

where we have applied Eq.(5). It is evident that gluon and electromagnetic penguin contri- 
butions to the CP-odd phase of the electric dipole DE amplitude are equally important, as 
it should be. 

At this point we would like to comment our work in relation to the study in Ref. [2] . Dib 
and Peccei first calculated CP asymmetry for charged kaon decay into two pions and then 
applied the CPT relation f] 

T(K + -> tt+tt ) + T(K + -> 7T + 7r°7) = Y{K~ -> t^tt ) + T{K~ -> 7r~7r°7) (26) 

to estimate Ap for K — > 7r ± 7r°7 decays. On the contrary, we compute CP- violating asymme- 
try directly for the radiative decays of K . Therefore, the strong-interaction phase difference 
necessary for generating CP-odd asymmetry is sin(5| — 5%) in our case, while it is sin(5 7 — Sq) 
in Ref. [2], where 5 7 is the strong-interaction phase shift for K ± — > tt^tt amplitudes involving 
a 7T7T7 intermediate state. Apart from this, it seems to us that the numerical discrepancy 
between the present work and Ref. [2] lies mainly in the fact that a factor of l/(47r) is missing 
in Eq.(15) of Ref. [2] for the effective Lagrangian of electromagnetic penguins. Consequently, 
Ar is overestimated by a factor of (47r) 2 ; in other words, the predicted upper bound for Ar 
in Ref. [2] should read 5.6 x 10~ 6 instead of 9 x 10~ 4 . Nevertheless, Dib and Peccei did point 
out the importance of the electromagnetic penguin diagram, which is no longer negligible for 
3 In principle, one should also include the decay rates of K — ► 3-7T decays, namely T(K + — * 7r + 7r + 7r~) + 
T{K + — » 7T + 7r 7r°), to the l.h.s. of Eq.(26) and their charge conjugate to the r.h.s. 



m t > Mw- From Eqs.(20) and (25) we see that gluon and electromagnetic penguin diagrams 
contribute constructively to 4>e- As a result, 

sin0£ = (sin0 £ ) g i uon + (sin0 E ) em = -1.4 x 10~ 4 sin<5i 3 . (27) 

It remains to work out the quantity z defined in Eq.(12). It can be recast in terms of the 
variables x = 2k ■ q/m 2 K and y = 2k ■ p+/m 2 K : 

z=\[x(l-y)-{l-y) 2 \. (28) 

Note that in the kaon rest frame, x = 2E 1 /rriK, y = 2E + /rriK- It is easily seen from Eq.(28) 
that the maximum z for a given photon energy E 1 in the cm. is given by 

(*)»« = 7 (— V- ( 29 ) 



4 \rriK 

Since (<5 1 - <Jg) ^ io° [15,16] and A(AT+ -> tt+tt ) = 1.829 x lO" 8 GeV [19], it follows from 
Eqs.(7), (12), (25) and (29) that 

0.75 x 1(T 5 ( mn E J, V Y 
A(E 7 ) = AmM^ , (30) 

! + 0-31 (too^v) 

This CP-odd asymmetry ranges from 2 x 1CT 6 to 1 x 1CT 5 when E 1 varies from 50 MeV to 
its highest value of 170 MeV. 

To conclude, we have shown in the 1/N C approach that the El amplitude of DE in 
— > 7r ± 7T°7 decays is not negligible. Therefore, a measurement of the interference of 
inner bremsstrahlung with electric dipole transitions is important for differentiating between 
various models. We also pointed out that CP- violating asymmetry between the Dalitz plots of 
— > 7r ± 7r°7 decays receive equally important contributions from gluon and electromganetic 
penguin diagrams. The magnitude of CP asymmetry ranges from 2 x 10 -6 to 1 x 10~ 5 when 
the cm. photon energy varies from 50 MeV to 170 MeV. 

Note added: After this work was typed, we learned a paper by G. Ecker, H. Neufeld 
and A. Pich [CERN-TH-6920/93 and UWThPh- 1993-22] in which the decay K+ -> 7r+7r 7 
is analyzed and a potentially sizeable electric amplitude interfering with bremsstrahlung is 
emphasized. 
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